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Abstract. The inelastic x-ray scattering spectrum for phonons of Ai-symmetry 
including the CuO bond-stretching phonon dispersion is analyzed by a Lorentz 
fit in HgBa2Cu04 and Bi2Sr2Cu06, respectively, using recently calculated phonon 
frequencies as input parameters. The resulting mode frequencies of the fit are almost 
all in good agreement with the calculated data. An exception is the second highest Ai- 
branch compromising the bond-stretching modes which disagrees in both compounds 
with the calculations. This branch unlike the calculations shows an anomalous 
softening with a minimum around the wavevector q = —(0.25, 0, 0). Such a disparity 
with the calculated results, that are based on the assumption of an undisturbed 
translation- and point group invariant electronic structure of the CuO plane, indicates 
some static charge inhomogenities in the measured probes. Most likely these will 
be charge stripes along the CuO bonds which have the strongest coupling to certain 
longitudinal bond-stretching modes that in turn selfconsistently induce corresponding 
dynamic charge inhomogenities. The symmetry breaking by the mix of dynamic and 
static charge inhomogenities can lead to a reconstruction of the Fermi surface into 
small pockets. 
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In recent work [U [2] the IXS phonon spectrum of optimally doped HgBa 2 Cu04 
and Bi 2 Sr 2 Cu0 6 , respectively, has been measured for the Ax-symmetry modes. The 
spectra have been interpreted by a Voigt fit function. For the interesting high 
frequency part of the spectra where the CuO bond-stretching modes (BSM) are 
expected that are known to display an anomalous softening upon doping in the cuprates 
jaUigiQHiaiaHniinra, onl y one Phonon peak has been fitted in case of HgBa 2 Cu0 4 
[1] and two in case of Bi 2 Sr 2 Cu0 6 [2]. 

On the other hand, from our microscopic calculation of the phonon dynamics of 
HgBa 2 Cu0 4 and Bi 2 Sr 2 Cu0 6 [13] we find in the high energy region of the IXS spectra 
along ^r-(C> 0, 0) up to four Ai modes as far as HgBa 2 Cu04 is concerned and up to five Ai 
modes in Bi 2 Sr 2 CuOg. These modes interact strongly with each other and constitute a 
highly nontrivial anticrossing scenario. As a consequence a multiple mode fit is necessary 
to obtain more reliable results. The strength of the measured peaks strongly decreases 
for larger £ values such that at C ~ 0.5 the peaks are very weak, see e.g. figure [I] 
and an unique interpretation taking only one or two peaks into account is incomplete 
because of the many modes of the same symmetry in this region of energy. Nevertheless, 
it is very interesting that in the limited Voigt fit of the spectra in [H [2] a hint to an 
anomalous softening of the BSM- dispersion is extracted at £ ~ 0.3 in the Hg compound 
and at £ ~ 0.25 in the Bi compound. Such an enhanced softening at these C values has 
not been detected in our calculations for the Ai modes. Only the presumably generic 
softening in the cuprates of the half-breathing mode (HBM) in HgBa 2 Cu04 and of two 
half-breathing modes in Bi 2 Sr 2 Cu06, generated by the strong nonlocal electron-phonon 
interaction (EPI) mediated by the charge fluctuations (CF) on the outer electron shells 
of the ions, has been found in our calculations [13] . 

In the presence of these theoretical results and the complex anticrossing of the Ax 
modes it is very desirable to perform a more complete analysis of the IXS measurements 
on the basis of the calculation of the dispersion of the Ai modes as given in [13] . In 
such a reconstruction of the IXS phonon spectra we do not apply the more complicated 
Voigt fit but use instead a more simple Lorentz fit function, which does the job very 
well. This can be concluded from figured] (a) where the IXS spectrum of HgBa 2 Cu04 
PQ is reconstructed for q = —(£,0,0) using as in pQ only one peak in the Lorentz fit. 
The phonon frequencies fitted at C = 0.11, 0.17, 0.23, 0.29, 0.36, 0.42 and 0.48 are 16.28 
(16.32), 15.76 (15.73), 14.61 (14.65), 13.38 (13.73), 13.15 (13.14), 13.40 (13.53) and 
13.51 (13.63) THz. The data within the brackets are the frequencies as obtained from 
the Voigt fit [1] . We recognize a very good agreement between both fitting procedures. 
As far as Bi 2 Sr 2 Cu06 is concerned we also get nearly undistinguishable results for a two 
mode fit in the two fitting schemes, so we take the more simple Lorentz fit in this work. 

In figure [1] (b) we present the results of the Lorentz fit for the analysis the IXS 
spectrum of HgBa 2 Cu04 within the framework of the calculated Ai-modes in the 
relevant energy range [13J. As already mentioned up to four mode frequencies must 
be considered. The calculated frequencies at the corresponding ( values are taken as 
input parameters of the Lorentz fit. From the results displayed in figured] (b) we find 
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Figure 1. (a) Reconstructed IXS phonon spectrum of HgBa2Cu04 from pQ for 
q = ^(£,0,0) using a Lorentz fit with one mode only (thick line), the corresponding 
Lorentz peaks of the modes are also shown (thin line) . The broken line represents the 
clastic background. The inelastic cross section N(E) is given in arbitrary units, (b) 
Same as in (a) but using a Lorentz fit on the basis of the calculated phonon dispersion 
for modes of Ai-symmetry in the energy range considered [15] . (c) Calculated phonon 
dispersion of HgBa2CuC>4 for modes of Ai-symmetry according to [T5J. The symbols 
■ and □ denote the phonon frequencies from the Lorentz fit and • the half-breathing 
mode. The vertical bars indicate the FWHM of the peaks determined in the fitting. 



that the IXS spectrum is well represented by the multiple mode fit. The Lorentz peaks 
as obtained by the fit are also shown in the figure for each ( value. 

The resulting phonon frequencies from the fit are shown in figure [TJ (c) together 
with the calculated dispersion of the A 1 modes [13]. With exception of the second 
highest Ax branch comprising the BSM that displays the anomalous softening starting 
at £ ~ | with a minimum around ( = |, the resulting mode frequencies are in reasonable 
agreement with the calculated branches. 

In case of Bi 2 Sr 2 Cu06 the Lorentz fit of the IXS spectrum on the basis of the 
calculated Ai-modes [13] is shown in figure [2] (b) and likewise as in HgBa2Cu04 a good 
resolution of the spectrum is achieved. The extracted phonon frequencies from the fit 
are shown in figure [2] (c). For a comparison the calculated dispersion of the Ax-modes 
[T3] are displayed in the figure. Leaving aside the second highest branch we find a 
good agreement of the calculated branches with the modes from the Lorentz fit. As for 
HgBa2Cu04 the second highest Ai branch displays an anomalous softening starting at 
C ~ g with a minimum around ( = \- 
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Figure 2. (a) Reconstructed IXS phonon spectrum of Bi2Sr2Cu06 from [2] for 
q = 0, 0) using a Lorentz fit with two modes only (thick line), the corresponding 
Lorentz peaks of the modes are also shown (thin line) . The broken line represents the 
clastic background. The inelastic cross section N(E) is given in arbitrary units, (b) 
Same as in (a) but using a Lorentz fit on the basis of the calculated phonon dispersion 
for modes of Ai-symmetry in the energy range considered [15] . (c) Calculated phonon 
dispersion of Bi2Sr2Cu6 for modes of Ai-symmetry according to [T5J. The symbols ■ 
and □ denote the phonon frequencies from the Lorentz fit and • the half-breathing 
modes. The vertical bars indicate the FWHM of the peaks determined in the fitting. 



Strong coupling of the BSM with dynamic charge stripes generated selfconsistently 
via strong nonlocal electron-phonon interaction (EPI) by the BSM themselves leading 
to an anomalous softening has already been found earlier in the cuprate based 
superconductors for the case of an undisturbed translation- and point group invariant 
electronic structure of the CuO plane without any charge inhomogenities [141 [T5l [161 
rm [181 [19] . The same holds true for HgBa 2 Cu0 4 and Bi 2 Sr 2 Cu0 6 as can be seen from 
the calculations in this work displayed in figure El where the charge redistribution 5p(r) 
for the half-breathing mode (HBM) (C = 5) induced by nonlocal EPI is shown. In the 
HBM propagating in A ~ (1, 0, 0) direction the O x ions move coherently along the CuO 
bond towards or away from the silent Cu ion. For the detailed displacement pattern, 
see [15] . 

Figure [3J illustrates that the HBM creates in both compounds dynamic charge 
ordering bp via charge fluctuations at the ions due to nonlocal EPI in form of localized 
stripes of alternating sign in the CuO plane. The period of the pattern is twice the 
lattice constant for £ = |. In case of Bi 2 Sr 2 CuOg these dynamic charge inhomogenities 
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Figure 3. Dynamical charge stripes Sp as induced by the half-breathing mode (with 
the lower frequency) in E^S^CuOg (a) and in HgBa2Cu04 (b). The arrows indicate 
the vibrating oxygen ions. Units are in 10~ 4 e 2 /ae- Regions of space where electrons 
are accumulated are characterized by full lines and those where electrons are pushed 
away by broken-dotted lines. 



are only shown for one of the two HBM modes found in the calculations namely the 
one with the lower frequency, see figure [2] (c), because the charge redistribution is very 
similar for the other HBM. For HgBa 2 Cu04 we obtain only one HBM, see figured] (c), 
which induces the dynamic stripe pattern in figure [3] (b). In the HBM the dynamic 
charge stripes point along the x- or y-axis, respectively, if the oxygen ions move along 
the CuO bond in y- or x-direction, respectively. 

It should be noted that for the oxygen breathing mode at the X point, see [13], the 
dynamic stripes excited point along the diagonals of the CuO plane. Furthermore, it 
should be remarked that dynamic charge inhomogenities like the dynamic stripe patterns 
would not be present in systems with only local EPI at work. For example in a high 
density homogeneous electron gas prevailing in simple metals. This is because of perfect 
screening of the changes of the Coulomb potential related to the displacement of the 
ions in a phonon mode. On the other hand, our calculations exhibit a strong nonlocal 
EPI in the cuprates essentially due to the poor screening in these materials with a large 
component of ionic binding, see. e.g. [19]. This fact demonstrates the special role of 
a strong nonlocal EPI for the physics of these compounds in their normal as well as in 
the superconducting state. 

Such a strong coupling is also underlined by recent experimental results for the 
quasi-particle dispersion in (Bi,Pb)2Sr 2 Cu06+<5 [20J. The super-high resolution laser- 
based angle-resolved photoemission spectroscopy measurements reported in [20] show 
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Figure 4. Calculated phonon density of states for Bi2Sr2Cu06- The gray domains 
mark the peaks of the Eliashberg spectral function [2U] . 



unambiguously that the ~ 70 meV nodal kink in the dispersion is due to strong EPI 
with multiple phonon modes as seen in the Eliashberg spectral function extracted from 
the measured real part of the electron self-energy. Remarkably, the peak structure of 
the Eliashberg function is in very good agreement with our calculated phonon density 
of states for Bi 2 Sr 2 Cu0 6 [13], see figure HI 

Quite generally, our calculations for the cuprates demonstrate the enhanced 
susceptibility of the BSM for inducing dynamic charge inhomogenities in form of 
stripes in the CuO plane. So, it can safely argued that the latter also will occur in 
a modified way in systems with static stripe charge order breaking translation- and 
point group invariance. The anomalous softening around C = | found by the analysis of 
the IXS phonon spectra then reflects the strong nonlocal coupling of the BSM with 
their selfconsistently induced dynamic charge inhomogenities (dynamic stripes) but 
now in probes with a static charge inhomogenity (static stripes) presumably present. 
Excess holes in a stripe will even enhance the strong nonlocal coupling to the BSM 
and consequently also the magnitudes of the CF's and the related mode softening is 
increased. An interrelation between the number of holes and softening is supported by 
the experiment and our calculations for the Cu-O-bond stretching modes [TH [15] . 

The symmetry breaking of the electronic structure by charge inhomogenity in form 
of static stripe patterns accompanied by phonon-induced dynamic fluctuating patterns 
shown by our calculations to occur very likely may lead to a corresponding reconstruction 
of the Fermi surface into small pockets. This reconstruction is in turn the origin for 
the quantum oscillations measured in some cuprates, see e.g. [2Tj. In this context it is 
interesting, that in [22] the quantum oscillations are attributed to symmetry breaking 
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by a certain type of stripe order in form of a combined charge/ spin modulation. Other 
authors attribute the small Fermi surface pockets to magnetic fluctuations [23], to d- 
density wave [21] or to an antiferromagnetic (n, 7r) spin density wave [25] . 

Finally, it should be added that in the literature dynamic charge inhomogenities 
related to electron-phonon coupling are discussed from the experimental side for 
Las-^Ba^CuC^, YBa 2 Cu 3 7 or Ba . 6 Ko.4Bi0 3 in [261 [2H1 [2H EQl EH [321 [33] and 
from the theoretical side e.g. in [31]. 

In summary, our analysis of the IXS phonon spectra for HgBa2Cu04 and 
Bi 2 Sr 2 Cu06 within the framework of calculated phonon dispersion curves points to the 
existence of dynamic charge inhomogenities generated by the BSM via strong nonlocal 
EPI of the CF-type. These modes excite dynamic charge order and consequently can 
also serve as indirect probes for corresponding static charge order in the cuprates. Thus, 
we have a mix of static and fluctuating charge order in the materials. 
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